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1. Introduction 



Pion polarizabilities are fundamental structure parameters characterizing the behavior of the 
pion in an external electromagnetic field. The dipole electric (cc\) and magnetic (Z^) pion polariz- 
abilities measure the response of the pion to quasistatic electric and magnetic fields. On the other 
hand, the quadrupole polarizabilities a 2 and /3 2 measure the electric and magnetic quadrupole 
moments induced in the pion in the presence of an applied field gradient. 

The generalized dipole (CC\ and j8i) and quadrupole (0:2 and /3 2 ) polarizabilities are defined 
[1,2] through expansion of the non-Born helicity amplitudes of Compton scattering on the pion in 
powers of t at fixed s = jx 2 



M ++ (s = IX 2 , t) = TtfX 

% 



2(a 1 -/3 1 ) + ^(a 2 -/3 2 ; 



M + _(s = i i I ,t) = - 2(a 1 +j3 1 ) + -(a 2 + /3 2 ) (1-1) 
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where s = (q\ + k]) 2 , t = (k\— fc 2 ) 2 , 12 and k\ , & 2 are the pion and photon four-momenta), and 
/X is the pion mass. In the following the dipole and quadrupole polarizabilities are given in units 
10~ 4 fm 3 and 10~ 4 fm 5 , respectively. 

It should be noted that the generalized pion polarizabilities very strongly differ from intrinsic 
polarizabilities a ( and j3 ; [3]. For example, the nonrelativistic expression of the generalized electric 
dipole polarizability a\ is equal to 

a\ = «i + ^jj <r l >=«i + 15.5, (1.2) 

where < > is the mean-square pion radius, a is the fine-structure constant. 

The expression for the generalized magnetic dipole polarizability j3i is more complicated. Ad- 
ditional contributions could be both paramagnetic and diamagnetic. So, it is difficult to determine 
the nature of magnetic susceptibility of j3i at present. 

In the following we will omit the word "generalized". 

The values of the pion polarizabilities are very sensitive to predictions of different theoretical 
models. Therefore, accurate experimental determination of these parameters is very important for 
testing the validity of such models. 

For example, the results of calculations of (ot\ — j8i)^± in Refs [2,4-6] are at variance with the 
predictions of chiral perturbation theory (ChPT) [7, 8]. On the other hand, P. Pasquini, D. Drechsel, 
and S. Scherer (PDS) in Ref. [9] claim that as the absorptive part of the Compton amplitudes in 
Refs [2, 4-6] is expressed by Breit-Wigner poles with energy dependent coupling constants and 
decay widths, there must appear additional spurious singularities. As a result, the values of the 
polarizabilities obtained in [2,4-6] have to be modified essentially. 

In the present paper we examine the statement of PDS and give an overview of the present 
situation in the field of investigation of the dipole and quadrupole polarizabilities of the charged and 
neutral pions in the frameworks of different DSRs and ChPT and compare the results of calculations 
with available experimental data. 
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2. Dispersion sum rules for dipole polarizabilities 

DSRs for the difference and sum of electric and magnetic pion polarizabilities have been con- 
structed using dispersion relations (DRs) for the helicity amplitudes M ++ and M + _, respectively. 
It has been shown in Ref. [10] that these amplitudes have no kinematical singularities or zeroes. 

In Refs [2,4, 6] DSR has been constructed for (cc\ — j3i) using DR for the amplitude M ++ at 
fixed u = jx 2 (where u = 2jx 2 — s — t). In this case, the Regge-pole model allows the use of DR 
without subtractions [10]. Such a DSR is 



f ImM ++ (t',u = a 2 ) dt' f ImM ++ (s'u = u 2 )ds'\ , n „ 

4/i 2 4/i 2 J 

The imaginary parts of the amplitudes in these DR and DSR are approximated by the contributions 
of meson resonances using Breit-Wigner expressions. For example, the contributions of the vector 
mesons are determined as 



for s > 4jU 2 and ImM^l(s,t) = for s < 4jx 2 , where 

3 




My x r v ^ y7l , r = f 5 4fi 4 ) M v r v . (2.3) 



My-H 2 / \M 2 -4n 



Here My, FV, and Ty^yn are the mass, full and yn decay widths of the vector mesons, respectively. 
A dependence of the width on the energy is conditioned by the threshold behavior. The energy 
dependence of the coupling constant gy appears via an expression for the total cross section of the 
process yn — ► yn through the vector meson contribution. 

In order to check the possibility of the appearance of additional singularities in our dispersion 
approach, we calculate the contributions of all mesons, except a, to our DSR by the zero-width 
approximation. 

The results of such calculations of (ai — /3i) z are listed in Tables 1, 2 together with the com- 
plete calculations of {a,\ - jSi)/ obtained in Ref. [2]. 



Table 1: The DSR predictions for («i - P\) % ±. 







P 


bi 


a\ 


ai 


fo 


fo 


a 


E 


AZ 




-h)f 


-1.15 


0.93 


2.26 


1.51 


0.58 


0.02 


9.45 


13.60 


2.15 


(0Cl 




-1.11 


0.85 


3.39 


1.51 


0.59 


0.03 


9.45 


13.70 





As seen from the Tables the zero-width approximation results practically coincide with the 
calculations of Ref. [2] which are beyond such an approximation. 

The coefficient 1/y^in the coupling constant of the a-meson amplitude in Ref. [2,4] provides 
the correct asymptotic behavior for the convergence of the integral over t in our DSR and does 
not lead to additional singularities. To check it we have calculated this integral using the energy 
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Table 2: The DSR predictions for {a x - fr)^. 



P (0 /o /o g 



E 



AE 



(a x — ft)/ -1.58 -12.56 -0.04 0.60 0.02 10.07 -3.49 2.13 
(ai-pi) z -1.99 -11.81 -0.04 0.61 O02 10.07 -3.14 



independent values of the decay width and the coupling constant of the a-meson. It has not lead 
to essential changes of the calculation results presented in the Tables. For example, the value of 
((Xi — p\) n ± z would be equal to 13.1. 

Besides, we compare our DSR calculation results in the ^-channel with the predictions of DSR 
obtained at the fixed angle Q yn = 180° [5]. In this case 



where v is the incident photon energy in the lab. system, a (yes) and a (no) stand for the sum of 
the photoabsorption cross sections containing, respectively, parity-flip and -nonflip multipoles. 

The best way to calculate this integral would be the use of experimental values for the cross 
sections. On the other hand, if these experimental data are well described by some function in the 
physical region of the process under consideration, then we can use this function, but only in this 
region, without continuation into unphysical regions. Usual Breit-Wigner forms for the photoab- 
sorption cross sections with energy dependent decay widths and correct asymptotic behavior are 
such functions. Therefore, we should use them only in the physical regions of these processes. If 
we add a contribution from spurious singularities, which are out of the physical region considered 
in DSR, we would have an additional contribution to the result obtained from the integration of 
the experimental cross section. This is a gross mistake. So, there are no problems with additional 
spurious singularities in the derivation and the calculation of Eq. (2.4). The results of the calcula- 
tion of this expression are very close to the values of ((X\ — jSi)/ given in Tables (1) and (2) for p, 
CO, <j), a\, and ai mesons which saturate the DSR integrals in the ^-channel. This result confirms 
the absence of additional singularities in our approach for the ^-channel integral of our DSRs. The 
^-channel contributions with / = / = are the same for both DSR at fixed u = jx 2 and DSR at fixed 
dy K = 180°. There are the same arguments why additional singularities are absent in the ^-channel 
too. 

It should be noted that in the work of S.S. Kamalov, L. Tiator, D. Drechsel et al. [11] the neutral 
pion photoproduction and electroproduction at the threshold were analyzed using DR. In this work 
the resonance contribution to the imaginary parts of the amplitudes of the pion photoproduction 
and electroproduction were given in terms of Breit-Wigner expressions with the energy dependent 
decay widths and coupling constants. According to the main statement of PDS, there are very many 
additional singularities in the resonance amplitudes considered in this work. The results of the 
calculations in the work [11] were obtained without consideration of any additional singularities 
and are in very good agreement with the experimental data for the pion photoproduction in the 
threshold region. However, if, according to PDS, one takes into account these singularities, this 
would lead to additional contributions and, as a result, to a disagreement with the experiment. This 
confirms that an account of such singularities is a mistake. 




3/1/2 



(2.4) 
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It is worth noting that the calculation of (a,\ — Pi)x ± in the framework of DSR at finite energy 
[12], which takes into account the ^-channel and the Regge-pole asymptotic contribution only, 
yielded (ot\ — j3i) ff ± = 10.3 ± 1.3. This value practically coincides with our result (see Table 1). 

As for the sum (ai + j3i), these values are calculated using Baldin's DSR 



(„ 4. fi ^ - 1 f a T(v)d v 



(2.5) 



if 



where Oj is the total cross section of the /^-interaction. These DSR results in 

(ai+jSi)«± = 0.166±0.024, 

(osi +j3i)«o =0.802 ±0.035. (2.6) 

On the other hand, two-loop ChPT calculations [7, 13] give 

(ai-ftV =5.7 ±1.0, 
(ai+ftV =0.16, 
(«i-j3 1 ) 7t o = -1.9±0.2, 

(ai+ft) a o = l.l±0.3. (2.7) 

So, the results of the ChPT calculations for the sum and difference of the dipole polarizabilities 
of 7T° and the sum for the charged pions do not conflict within the errors with predictions of DSRs. 

Let us consider possible reasons of the discrepancy between the predictions of DSRs and ChPT 
for (oci — j3i)jr±. The main contribution to the DSRs for (ai — pi) n ± is given by the a-meson. 
However, this meson is taken into account only partially in the present ChPT calculations. 

Consider the methods of the calculation of the vector meson contribution in the frameworks of 
DSRs and ChPT. In the narrow width approximation we have from Eq. (2.2) 

ImM^l(s,t) = -- g 2 v s8{s-Ml). 
7T 



Then the DSR calculation gives 

•oM , , ( v — a 2 t = n\ = 

(Ml -v 2 



ReM ++ ( s = ii 2 ,t = 0)= , a / v v ? _ (2. 8) 



In the case of ChPT the authors of Ref. [7] used 

Re « ++ ( S = ,,V = 0) = ^!_. (2.9) 

The absolute value of the amplitude (2.9) is smaller than (2.8) by a factor My /fx 2 . From the point 
of view of analyticity, the result (2.9) could be obtained if DR with one subtraction at s = is used 
for the amplitude M ++ (s,t). However, an additional subtraction constant M ++ (s = 0,t = 0) then 
appears, which was not considered in the available ChPT calculations. 

In the case of the difference of the dipole polarizabilities of the 7T°-meson, the big contribution 
of the a-meson to DSR is cancelled by the big contribution of the co-meson. On the other hand, in 
the ChPT calculations the a-meson is only partially included and the ft)-meson gives a very small 
contribution to this difference. As a result, the DSR and ChPT predictions for (cc\ — /3i) w o are rather 
close. 
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Table 3: The experimental data presently available for (<X\ —j5i) K ±. In [16-20] («i — was determined 
by using the constraint a ljc ± = — /3 1?I ±. 



Experiments 


(ai-ft),^ 


yp -► JK+n MAMI (2005) [14] 


1 1.6 ± l.5 stat ± 3-0 S y S t ± 0.5 mo ^ 


yp — »■ 77T + « Lebedev Phys. Inst. (1984) [15] 


40 ±20 


7TA -»• 77r"A Serpukhov (1983) [16] 


13.6±2.8±2.4 


tit A -► yjrA COMPASS (2007) [17] 


5 .0 ± 3 .4 (preliminary) 


77 ^ 71+ 71" (£ y < 700 MeV) 




D. Babusci etal. (1992) [18] 




PLUTO [21] 


38.2±9.6±11.4 


DM 1 [22] 


34.4 ±9.2 


DM 2 [23] 


52.6 ±14.8 


MARK II [24] 


4.4 ±3.2 


J.F. Donoghue, B.R. Holstein (1993) [19] 


5.4 


MARK II [24] 




A.E. Kaloshin, V.V. Serebryakov (1994) [20] 


5.25 ±0.95 


MARK II [24] 




L.V. Fil'kov, V.L. Kashevarov (2006) [4] 


n +2.6 


77 -> 7T+ 7T~ fit of data [24-29] 




from threshold to 2.5 GeV 





3. Experimental data for dipole polarizabilities of charged pions 

By now the values of the pion polarizabilities were determined by analyzing the processes 
7t~A — > 77T~A, 7/7 — > yn + n, and 77 — ► 7T7T. The experimental information available so far for the 
difference of the dipole polarizabilities of charged pions is summarized in Table 3. 

The values of the experimental cross sections of the process 77^ K + n~ in the energy region 
Ey < 700 MeV are very ambiguous. As a result, the values of (cc\ — j3i) w ±, obtained from analyses 
of these data, lie in the interval 4.4-52.6. The analyses of the data of Mark II [24] only have given 
(CC\ — pi) n ± close to the ChPT result. 

The difference [a,\ — j3i) w ± found from the global fit to all available experimental data of the 
process 77 — > % + %~ in the energy region from the threshold to 2500 MeV [2] agrees very well 
with the results [16] obtained from the scattering of high energy %~ mesons off the Coulomb field 
of heavy nuclei and from the radiative photoproduction of n + from the proton at MAMI [14] and 
in Lebedev Physical Institute [15] (see Table 3) and with the DSR calculations. However, these 
values of (<X\ — /3i) w ± deviate essentially from the ChPT calculations [7, 8]. 

Results of polarizability extraction from the process n~A — > yn~A strongly depend on the 
momentum transfer Q 2 . In this reaction the Coulomb amplitude dominates for Q 2 < 10~ 4 (GeV/c) 2 . 
In the region of Q 2 ~ 10~ 3 (GeV/c) 2 Coulomb and nuclear contributions are of similar size. In this 
region the nuclear contribution, in particular an interference between the Coulomb and nuclear 
amplitudes, should be taken into account. In the work [16] the authors considered Q 2 < 6 x 10~ 4 
(Gev/c) 2 , while the authors of Ref. [17] worked at Q 2 < 7.5 x 10~ 3 (Gev/c) 2 . In this region of Q 2 the 
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Table 4: The quadrupole polarizabilities of the neutral and charged pions. 









ChPT [7,13] 




fit [2,4] 


DSR [2] 


to one-loop 


to two-loops 




39.7 ±0.02 


39.72 ±8.01 




37.6 ±3.3 




-0.181 ±0.004 


-0.171 ±0.067 




0.04 


(afc-ft)** 


25.01JB 


25.75 ±7.03 


11.9 


16.2 [21.6] 


(a 2 + j3 2 ) ff± 


0.133 ±0.015 


0.121 ±0.064 





-0.001 [-0.001] 



contribution of the interference between the Coulomb and nuclear amplitudes is very large [30,31]. 
However, it was not taken into account in the work [17]. This is the main reason of the difference 
between the Serpukhov and COMPASS results. Moreover, for the total energy, in the y% c.m.s., 
W > 450 MeV and d r y ~ 180°, the a-meson contribution should be taken into account [6]. 

4. Pion quadrupole polarizabilities 

DSRs for the difference and the sum of the quadrupole polarizabilities have been obtained 
with the help of DRs at fixed u = ji 2 with one subtraction for the amplitudes M++ and M+_, 
respectively: 



(«2-/3 2 ) = - f - 



oo oo 

f ImM ++ (t',u = LI 2 ) dt' r ImM ++ (s',u = LI 2 ) ds' 

.1 T 2 J (s'-v 2 ) 2 

4/I 2 4/i 2 



(4.1) 



(a 2 + /3 2 ) = ^ 



ImM + -(t',u = Ll) dt 1 f ImM + ^(s',u = Ll ) ds' 



t' 2 J {s'-ll 2 ) 2 

4/i 2 4/i 2 



_ r ImM+_(s',u = Ll^ds' I 
J is' LI 2 ) 2 | 



The corresponding values of the quadrupole polarizabilities have been found in Refs [2, 4] 
by fitting the experimental total cross sections of the processes 77 — > n°n° [32, 33] and 77 — > 
7T + 7T~ [24-29] in the energy regions from the thresholds to 2.25 and 2.5 GeV, respectively. The fit 
functions were constructed by using DRs with subtractions for the amplitudes M ++ and M + _. 

The values of the quadrupole polarizabilities found in Refs. [2,4] and the predictions of DSRs 
[2] and ChPT [7, 13] are listed in Table 4. The numbers in brackets correspond to the order p 6 low 
energy constants from Ref. [34]. As seen from this Table, all values of the polarizabilities found in 
Ref. [2,4] are in good agreement with the DSR predictions [2]. 

The difference of the quadrupole polarizabilities (a 2 — /fe)^ obtained in Refs. [2,4] disagrees 
with the present two-loop ChPT calculations [7, 13]. One of the sources of this disagreement is 
the bad knowledge of the low energy constants. Moreover, it should be noted that in this case 
the two-loop contribution generates nearly 100% as compared to the one-loop result. The ChPT 
calculations of (a 2 + p 2 ) give an opposite sign. However, calculations of (a 2 + j3 2 ) at order p 6 
determine only the leading order term in ChPT. Therefore, contributions at p 8 could be essential, 
and considerably more work is required to put the chiral prediction on a firm basis in this case [7]. 

It is worth noting that calculations of the dipole and quadrupole pion polarizabilities in the 
frame of the Nambu-Jona-Lasino model [35] agree within errors with the DSR [2,4] predictions. 
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5. Summary 

We showed that there are no problems with additional spurious singularities in the DSRs and 
DRs considered. The difference between the predictions of the DSRs and ChPT for (oc\ — pi) n ± 
remains. This discrepancy is connected with a different account of the contribution of the o and 
vector mesons in the DSR and ChPT calculations. The disagreement between the DSR and ChPT 
predictions of the quadrupole polarizabilities is connected, in particular, with the bad knowledge 
of the low energy constants. Substantial corrections to the values of the quadrupole polarizabilities 
are expected from three-loop calculations. 

The authors thank A.I. L'vov and V.A. Petrun'kin for useful discussions. This research was 
supported by the DFG-RFBR (Grant No. 09-02-91330). 
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